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Abstract. 

LDA+DMFT (Local Density Approximation combined with Dynamical Mean- 
Field Theory) computation scheme has been used to calculate spectral properties 
of LaFeAsO - the parent compound for new high-Tc iron oxypnictides. Coulomb 
repulsion U and Hund's exchange J parameters for iron 3d electrons were calculated 
using first principles constrained density functional theory scheme in Wannier functions 
formalism. Resulting values strongly depend on the number of states taken into account 
in calculations: when full set of 0-2p, As-4p, and Fe-3(i orbitals with corresponding 
bands are included, computation results in U =3^4 eV and J=0.8 eV. In contrast to 
that when the basis set is restricted to Fe-3(i orbitals and bands only, computation gives 
much smaller parameter values F°=Q.8 eV, J=0.5 eV. However, DMFT calculations 
with both parameter sets and corresponding to them choice of basis functions result 
in weakly correlated electronic structure that is in agreement with experimental X-ray 
and photoemission spectra. 
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Figure 1. Total and partial densities of states for LaFeAsO obtained in DFT 
calculation in frame of LMTO method. 

Recent discovery of high-Tc superconductivity in iron oxypnictides LaOi_a;Fa;FeAs |1] 
has stimulated an intense experimental and theoretical activity. In striking similarity 
with high-Tc cuprates, undoped material LaFeAsO is not superconducting with antifer- 
romagnetic commensurate spin density wave developing below 150 K [2]. Only when 
electrons (or holes) are added to the system via doping, antiferromagnetism is suppressed 
and superconductivity appears. As it is generally accepted that Coulomb correlations 
between copper 3d electrons are responsible for cuprates anomalous properties, it is 
tempting to suggest that the same is true for iron 3d electrons in LaFeAsO. 

The ratio of Coulomb interaction parameter U and band width W determines 
correlation strength. If U/W < 1 then the system is weakly correlated and results 
of the Density Functional Theory (DFT) calculations are reliable enough to explain 
its electronic and magnetic properties. However, if U value is comparable with W or 
even larger then the system is in intermediate or strongly correlated regime and Coulomb 
interactions must be explicitly treated in electronic structure calculations. The partially 
filled bands formed by Fe-3d states in LaFeAsO have the width f«4 eV (see shaded area 
in the lower panel of Fig. [1]), so the estimation for Coulomb interaction parameter U 
should be compared with this value. 

In practical calculations, U is often considered as a free parameter to achieve 
the best agreement of calculated and measured properties of investigated system. 
Alternatively, U value could be estimated from the experimental spectra. However, 
the most attractive approach is to determine Coulomb interaction parameter U value 
in first principles non-empirical way. There are two such methods: constrained 
DFT scheme |3l H], where li-orbital occupancies in DFT calculations are fixed to the 
certain values and U is numerically determined as a derivative of d-orbital energy 
over its occupancy, and Random Phase Approximation (RPA) method [5], where 
screened Coulomb interaction between (i-electrons is calculated via perturbation theory. 
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Figure 2. (Colour online) Module square of drc2_y2-\i'ke Wannier function computed 
for Fe-3(i bands only (left panel) and for full set of 0-2p, As-4p and Fe-Sd bands (right 
panel). Big sphere in the center marks Fe ion position and four small spheres around 
it correspond to As neighbors. 

Recently, the RPA calculations for Coulomb interaction parameter U in LaFeAsO were 
reported [6], where U value was estimated as 1.8^2.7 eV. In Ref. [?] it was proposed to 
use in LaFeAsO ?7=4 eV obtained in RPA calculations for metallic iron [8]. 

This value for Coulomb parameter (with Hund's exchange parameter J =0.7 eV) was 
used in Dynamical Mean-Field Theory (DMFT) [9] calculations for LaFeAsO [Tl lTTlfTO] . 
Results of these works show iron 3d electrons being in intermediate or strongly correlated 
regime, as it is natural to be expected for Coulomb parameter value f/=4 eV and Fe-Sd 
band width ^4 eV. 

The most direct way to estimate correlation effects strength in a system under 
consideration is to compare the experimental spectra with densities of states (DOS) 
obtained in DFT calculations. For strongly correlated materials additional features in 
the experimental photoemission and absorption spectra appear that are absent in DFT 
DOS. Those features are interpreted as lower and upper Hubbard bands. If no Hubbard 
bands are observed and DOS obtained in DFT calculations satisfactorily describes the 
experimental spectra then the material is in weakly correlated regime. 

LaFeAsO was studied by soft X-ray absorption and emission spectroscopy fl2\ . 
X-ray absorption (O i^-edge) spectroscopy [13], and photoemission spectroscopy [H]. 
In all these works the conclusion was that DOS obtained in DFT calculations 
gave good agreement with the experimental spectra and the estimations |13] for 
Coulomb parameter value are U <1 eV. Such contradiction with results of DMFT 
calculations [7[ [HI HO] using U=4 eV shows that first principles calculation of Coulomb 
interaction parameter U value for LaFeAsO is needed to determine the correlation effects 
strength in this material. Results of such calculations by constrained DFT calculations 
are reported in the present work. We have obtained the value U <1 eV for calculation 
with basis restricted to Fe-3d bands in agreement with the estimates from spectroscopy. 

It is important to note that Coulomb interaction parameter U value depends on 
the choice of the model and, more specifically, on the choice of the orbital set that is 
used in the model. For example, in constrained DFT calculations for high-Tc cuprates 
the resulting U value for Cu ci-shell was found [I5] between 8 and 10 eV. The U value 
in this range was used in cluster calculations where all Cu d-orbitals and p-orbitals of 
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neighboring oxygens were taken into account and calculated spectra agree well with 
experimental data \1Q\. However, in one band model, where only — y"^ orbital per 
copper atom is explicitly included in the calculations [IT], the U value giving good 
agreement with experimental data falls down to 2.5-i-3.6 eV, that is 3-4 times smaller 
than constrained DFT value. 

The same situation occurs for titanium and vanadium oxides: the U value from 
constrained DFT calculations is ~6 eV and cluster calculations where all (i-orbitals 
and p-orbitals of neighboring oxygens were taken into account with U close to this 
value gave good agreement between calculated and experimental spectra [18]. However, 
in the model where only partially filled t2g orbitals are included, much smaller U 
value (corresponding to Slater integral F°=3.5 eV) gives the results in agreement with 
experimental data [l9] . 

It is interesting that such a small U value can be obtained in constrained DFT 
calculations for titanates and vanadates where only t2g-orbital occupancies are fixed 
while all other states (cg-orbitals of vanadium and p-orbitals of oxygens) are allowed to 
relax in self-consistent iterations [20l [T9] . So the calculation scheme used in constrained 
DFT (the set of the orbitals with fixed occupancies) should be consistent with basis set 
of the model where the calculated U value will be used. 

Another source of uncertainty in constrained DFT calculation scheme is a definition 
of atomic orbitals whose occupancies are fixed and energies are calculated. In some DFT 
methods, like Linearized Muffin-Tin Orbitals (LMTO), these orbitals could be identified 
with LMTO. However, in other DFT calculation schemes, where plane waves are used 
as a basis, like in pseudopotential method one should use more general definition for 
localized atomic like orbitals such as Wannier functions [23] (WFs). The practical way 
to calculate WFs for specific materials using projection of atomic orbitals on Bloch 
functions was developed in Ref. [24| . 

In Fig. [1] the total and partial DOS for LaFeAsO obtained in LMTO calculations 
are shown. Crystal field splitting for Fe-3d orbitals in this material is rather weak 
(Ac/=0.25 eV) and all five d orbitals of iron form common band in the energy region 
(—2, +2) eV relative to the Fermi level (see grey region on the bottom panel in Fig. [1]). 
There is a strong hybridization of iron t2g orbitals with p orbitals of arsenic atoms which 
form nearest neighbors tetrahedron around iron ion. This effect becomes apparent in 
the energy interval (—3, —2) eV (white region on the bottom panel in Fig. [1]) where 
band formed by p orbitals of arsenic is situated. More week hybridization with oxygen 

Table 1. The constrained DFT calculated values of average Coulomb interaction U 
and Hund's exchange J (eV) parameters for d-symmetry Wannier functions computed 
with two different sets of bands and orbitals. 



DFT method 


restricted to Fe-Sd bands 


full bands set 


TB-LMTO-ASA 
PWSCF 


f7=0.49, J=0.51 
[7=0.59, J=0.53 


f7=3.10, J=0.81 
f7=4.00, J=L02 
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Figure 3. Partial densities of states for Fe-3d orbitals obtained within the DPT (filled 
areas) and LDA+DMPT orbitally resolved spectral functions for "restricted basis" and 
F°=0.8 eV, J=0.5 eV (bold lines). 



p states reveals in (—5.5, —3) eV energy window (black region on the bottom panel in 
Fig.E). 

We have calculated Coulomb interaction U and Hund's exchange J parameters 
for WFs basis set via constrained DFT procedure with fixed occupancies for WFs of 
d symmetry. For this purpose we have used two calculation schemes based on two 
different DFT methods. One of them involves linearized muffin-tin orbitals produced by 
the TB-LMTO-ASA code [2Tj; corresponding WFs calculation procedure is described in 
details in Ref . [25j . The second one is based on the pseudopotential plane- wave method 
PWSCF, as implemented in the Quantum ESPRESSO package [22] and is described 
in Ref. [26]. The difference between the results of these two schemes could give an 
estimation for the error of U and J determination. 

The WFs are defined by the choice of Bloch functions Hilbert space and by a set 
of trial localized orbitals that will be projected on these Bloch functions [25]. We 
performed calculations for two different choices of Bloch functions and atomic orbitals. 
One of them includes only bands predominantly formed by Fe-3d orbitals in the energy 
window (—2, +2) eV and equal number of Fe-3d orbitals to be projected on the Bloch 
functions for these bands. That choice corresponds to the model where only five (i-orbital 
per Fe site are included but all arsenic and oxygen p-orbitals are omitted. Second choice 
includes all bands in energy window (—5.5, +2) eV that are formed by 0-2p, As-Ap and 
Fe-3d states and correspondingly full set of 0-2p, As-Ap and Fe-3d atomic orbitals to be 
projected on Bloch functions for these bands. That would correspond to the extended 
model where in addition to (i-orbitals all p-orbitals are included too. 

In both cases we obtained Hamiltonian in WF basis that reproduces exactly bands 
formed by Fe-3d states in the energy window (—2, +2) eV (Fig. [1]), but in the second 
case in addition to that bands formed by p-orbitals in the energy window (—5.5, —2) eV 
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Figure 4. Partial densities of states for Fe-Sd , As-4p and 0-2p states obtained within 
the DFT (fiUed areas) and LDA+DMFT orbitally resolved spectral functions for "full 
basis" calculations with eV, J=0.81 eV (bold lines). 



will be reproduced too. However, WFs with d-orbital symmetry computed in those two 
cases have very different spatial distribution. In Fig. [2] the module square of dx2_y2- 
like WF is plotted. While for the case when full set of bands and atomic orbitals was 
used (right panel) WF is nearly pure atomic d-orbital (iron states contribute 99%), WF 
computed using Fe-3d bands only is much more extended in space (left panel). It has 
significant weight on neighboring As ions with only 67% contribution from central iron 
atom. 

The physical reason for such effect is p-d hybridization that is treated explicitly in 
the case where both p- and (i-orbitals are included. In the case where only Fe-3d bands 
are included in calculation p-d hybridization reveals itself in the shape of WF. Fe-3d 
bands in the energy window (—2, +2) eV correspond to antibonding combination of 
Fe-3d and As-4p states and that is clearly seen on the left panel of Fig. [2l 

The different spatial distribution for two WFs calculated with full and restricted 
orbital bases can be expected to lead to different effective Coulomb interaction for 
electrons occupying these states. The results of constrained DFT calculations of the 
average Coulomb interaction U and Hund's exchange J parameters for electrons on 
WFs computed with two different set of bands and orbitals (and using two different 
DFT methods: LMTO and pseudopotential) are presented in Tab. [H 

One can see that very different Coulomb interaction strength is obtained for 
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Figure 5. (Colour online) Real (upper panel) and imaginary (lower panel) parts of 
LDA+DMFT self energy interpolated on real axis with the use of Fade approximant 
for "restricted basis" and F°^0.8 eV, J=0.5 eV. 

separate Fe-3d band and full bands set calculations. While the latter gives value 3-=-4 eV, 
close to previously used values [TJ [TTl [10], Fe-3d band restricted calculation results in 
0.5-i-0.6 eV, that is much smaller but agrees with spectroscopy estimations [13]. 

The main reason for such a drastic difference between two calculations results is 
very different spatial extension of the two WFs (see Fig. [2]): nearly complete localization 
on central iron atom for "full bands set" WF (99%) and only 67% for "Fe-Sc? band set" 
WF with 33% of WF on neighboring As atoms. Another reason for strong reduction of 
calculated If value in going from "full bands set" to "Fe-3(i band set" WF is additional 
screening via p-d hybridization with As-4p band that is situated just below Fe-3d band 
(see Fig. [1]). The effect of decreasing of the effective U value in several times going 
from full orbital model to restricted basis was found previously for high-Tc cuprates 
(f/=8-i-10 eV for full p-d-orbitals basis [15] and 2.5-^3.6 eV for one-band model [T7]). 

In constrained DFT calculations one obtains an average Coulomb interaction U 
that can be estimated[l] as f7 = F'' — J/2 . Hence, Slater integral F° can be calculated 
as FO = U + J/2. For "Fe-3rf band set" WF that gives F°=0.8 eV at J=0.5 eV. Coulomb 
parameters for "full basis" were calculated in the same way using data from the first row 
in Tab. [H, and were taken as F°=3.5 eV and J=0.81 eV. With this set of parameters 
we performed the LDA+DMFT [27] calculations (for detailed description of the present 
computation scheme see Ref. [25]). The DFT band structure was calculated within the 
TB-LMTO-ASA method [2T]. Crystal structure parameters were taken from Ref. [I]. 

The LDA+DMFT calculations were performed for both models: with restricted to 
Fe-3d states basis and with full basis including also As-p and 0-p states. For the later 
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Figure 6. (Colour online) Real (upper panel) and imaginary (lower panel) parts of 
LDA+DMFT self energy interpolated on real axis with the use of Fade approximant 
for "fuU basis" and F°^3.5 eV, J=0.81 eV. 

case double counting term U{nDMFT — ^) was used to obtain noninteracting Hamiltonian 
[26] . Here hdmft is the total number of (i-electrons obtained selfconsistently within the 
LDA+DMFT scheme. The effective impurity model for the DMFT was solved by the 
QMC method in Hirsh-Fye algorithm [28]. Calculations for both cases were performed 
for the value of inverse temperature (3=10 eV~^. Inverse temperature interval < t < (3 
was divided into 100 slices. 6 ■ 10® QMC sweeps were used in self-consistency loop within 
the LDA+DMFT scheme and 12 ■ 10® of QMC sweeps were used to calculate the spectral 
functions. 

The iron 3d orbitally resolved spectral functions obtained within DFT and 
LDA+DMFT calculations for "restricted basis" with ^"^=0.8 eV, J=0.5 eV are presented 
in Fig. [3l The influence of correlation effects on the electronic structure of LaFeAsO 
is minimal: there are relatively small changes of peak positions for 3z^ — r^, xy and 
— y"^ orbitals (the shift toward the Fermi energy) and practically unchanged picture 
of spectral function distribution for yz, zx bands. There is no appearance of either 
Kondo resonance peak on the Fermi level or Hubbard bands in the energy spectrum 
with such small values of U and J. 

Results for "full basis" LDA+DMFT calculations are presented in Fig. HI Note, 
that though much larger Coulomb U and J parameters (^"=3.5 eV, J=0.81 eV) were 
used in this calculation the spectra around the Fermi energy are very similar to those 
obtained in "restricted basis" calculations (Fig. [3]). The general shape of spectra does 
not show either Kondo resonance peak on the Fermi level or Hubbard bands and the 
features in Fe-d spectral functions below -2 eV correspond to hybridization with As-p 
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Figure 7. (Colour online) Splitting of Fe-d orbitals obtained in LDA and LDA+DMFT 
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Figure 8. (Colour online) Calculated Fe-d LDA+DMFT ("full basis" and F°=3.5 eV, 
J=0.81 eV) spectral function (green) and experimental Fe L3 XES spectrum (black 
circles) from Ref. |12) . 



and 0-p bands. The reason for such weak correlation effects in spite of relatively large 
Coulomb repulsion parameters is very strong hybridization of Fe-d orbitals with As- 
p states (see peaks in Fe-d spectral function in the region -2+-3 eV corresponding to 
admixture to As-p bands). This hybridization results in additional very effective channel 
of screening for Coulomb interaction between Fe-d electrons. 

This agrees with the results of soft X-ray absorption and emission spectroscopy 
study ^12]. It was concluded there that LaFeAsO does not represent strongly correlated 
system since Fe L3 X-ray emission spectra do not show any features that would indicate 
the presence of the low Hubbard band or the quasiparticle peak that were predicted by 
the LDA-I-DMFT analysis [71 [HI [10] with the large 11=4: eV. A comparison of the X-ray 
absorption spectra (O -ft"-edge) with the LDA calculations gave [13] an upper limit of the 
on-site Hubbard U ^1 eV. Photoemission spectroscopy study of LaFeAsO suggests [H] 
that the line shapes of Fe 2p core-level spectra correspond to an itinerant character of 
Fe 3d electrons. It was demonstrated there that the valence-band spectra are generally 
consistent with band-structure calculations except for the shifts of Fe Sd-derived peaks 
toward the Fermi level. Such a shift is indeed observed in our LDA-I-DMFT results 
(Fig. [3]). 
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Figure 9. (Colour online) Calculated total LDA+DMFT ("full basis" and F°=3.5 eV, 
J=0.81 eV) spectral function (green) and experimental NdFcAs02^Fi_2, PES spectrum 
(black circles) from Ref. [30] . 

The behavior of real part of self energy near zero frequency S(lo')|i^^o could give an 
important information about band narrowing and renormalization of electron mass in 
the system under consideration. Pade approximant [29] was used to obtain self energy 
on real frequencies. Results for calculation with both bases are presented in Figs. [Hi 
The calculated value of quasiparticle renormalization amplitude Z = {1 — |^^o))~^ 
was found to be 0.77, 0.63, 0.71, 0.46 for "restricted basis" and 0.56, 0.54, 0.45, 
0.56 for "full basis" for d^y, dyz (or dzx), d^^^^r'^, d^^-y^ orbitals, respectively. The 
values for "restricted basis" agree well with the effective narrowing of the LDA+DMFT 
spectral functions comparing with LDA DOS curves (Fig. [3]). Note, that quasiparticle 
renormalization amplitude Z = {1 — ^^^\^^q^)~^ from full basis calculations can not 
be directly used to estimate effective narrowing of the bands. Due to the explicit 
hybridization with As-p and 0-p states actual narrowing should be weaker. Keeping 
in mind the above mentioned, we can conclude that both models agree well with each 
other. The effective masses m* = Z^^ are 1.3, 1.59, 1.41, 2.17 ("restricted basis") for 
dxy, dyz (or dzx), d^z^-r^, dx2_y2 orbitals, respectively. The d^2_y2 orbital has the largest 
effective mass and exhibits most evident narrowing of LDA spectrum (see Fig. [3]). This 
orbital has its lobes directed to the empty space between nearest iron neighbors in iron 
plane. Hence it has the weakest overlapping, the smallest band width, and the largest 
U/W ratio. 

Small effective mass enhancement value shows that LaOFeAs belongs to the weakly 
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correlated systems in contrast to results of the LDA+DMFT calculations [7| by Haule 
et al where a strongly renormalized low energy band with a fraction of the original 
width {Z ^ 0.2-0.3) was found while most of the spectral weight is transferred into 
a broad Hubbard band at the binding energy ^4 eV. Authors of Ref. [7j report that 
"slightly enhanced Coulomb repulsion {U= 4.5 eV) opens the gap" so that the system 
is in strongly correlated regime on the edge of metal-insulator transition. 

In LaFeAsO the iron ion is in tetrahedral coordination of four As ions with slight 
tetragonal distortion of the tetrahedron. In tetrahedral symmetry group five d- 
orbitals should be split by crystal field [j| on low-energy doubly degenerate set 3z'^ — r^, 
xy corresponding to irreducible representation Cg and high-energy triply degenerate set 
— y^, xz, yz for representation t2g- We have calculated Wannier functions energy and 
have found that t2g-eg crystal field splitting parameter is very small Ac/ ~0.25 eV. The 
slight tetragonal distortion of the tetrahedron gives additional splitting of t2g and Cg 
levels with the following values for orbital energies (energy of the lowest Sz"^ — orbital 
is taken as a zero): e322_^2=0.00 eV, exy=0.03 eV, exz,yz=0.26 eV, £:^.2_y2=0.41 eV. 
Correlation effects lead not only to narrowing of the bands but also to substantial shift 
of Fe-d orbitals energies. For "full basis" calculation adding of i?e(E(0)) to the LDA 
orbital energies results in e3z2-r^=0.00 eV, exy=-0.37 eV, exz,yz=0-^0 eV, ex2_y2=0.20 eV 
(see Fig. [7]). Note that the actual splitting should be smaller due to p-d hybridization. 

Comparisons of spectral function calculated within LDA-I-DMFT for full basis 
set with various experimental spectra are presented in Figs. [H [91 One can see a 
good agreement in comparison of calculated Fe-d spectral function with Fe L3 XES 
spectrum (black circles) from Ref. [12]. The shoulder in experimental curve near -2.5 
eV correspond to the low energy peak in calculated spectrum that appears due to strong 
hybridization between Fe-d and As-p states (see also Fig. [1]). Note that this features in 
calculated curve is absent in the calculations with restricted basis (see orbital resolved 
spectra in Fig. [3]). 

The total spectrum (Fig. [9]) is in reasonable agreement with experimental 
photoemission data. The sharp peak at the Fermi energy corresponds to partially filled 
Fe-d band. Broad bands between -2 and 6 eV represent oxygen and arsenicum p bands. 

In conclusion, we have calculated the values of Coulomb interaction parameters 
U and J via constrained DFT procedure in the basis of Wannier functions. For 
minimal model including only Fe-3c? orbitals and bands we have obtained Coulomb 
parameters F°=0.8 eV, J=0.5 eV. For "full basis" calculations Coulomb parameters are 
F0=3.5 eV, J=0.8 eV. The LDA+DMFT calculation for both models with calculated 
parameters results in weakly correlated nature of iron d bands in this compound. This 
conclusion is supported by spectroscopic investigations of this material. As the choice 
of the orbitals to construct the effective Hubbard model from LDA results is to some 
extent arbitrary, showing that different choices leading to largely different interaction 
parameters yield similar physical observables provides an important consistency test of 

f The coordinate axes x,y in LaOFeAs crystal structure are rotated on 45 degrees from the standard 
tetrahedral notation so that xy and — orbitals are interchanged. 
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the whole procedure. 

Support by the Russian Foundation for Basic Researcli under Grant No. RFFI-07- 
02-00041, Civil Research and Development Foundation together with Russian Ministry 
of science and education through program Y4-P-05-15, Russian president grant for young 
scientists MK-1184.2007.2, President of Russian Federation fund of support for scientific 
schools grant 1941.2008.2 and Dynasty Foundation is gratefully acknowledged. J.K. 
acknowledges the support of SFB 484 of the Deutsche Forschungsgemeinschaft. 



References 



Kamihara Y, Watanabe T, Hirano M and Hosono H 2008 J. Am. Chem. Soc. 130 3296 

de la Cruz C, Huang Q, Lynn J W, Li J, Ratcliff II W, Zarestky J L, Mook H A, Chen G F, Luo 

J L, Wang N L and Dai P 2008 Nature 453 899 
Dederichs P H, Bliigel S, Zeller R and Akai H 1984 Phys. Rev. Lett. 53 2512 Gunnarsson O, 

Andersen O K, Jepsen O and Zaanen J 1989 Phys. Rev. B 39 1708 
Anisimov V I and Gunnarsson O 1991 Phys. Rev. B 43 7570 

Solovyev I V and Imada M 2005 Phys. Rev. B 71 045103 Aryasetiawan F, Karlsson K, Jepsen O, 

and Schonberger U 2006 Phys. Rev. B 74 125106 
Nakamura K, Arita R and Imada M, J. Phys. Soc. Japan 77 093711 (2008). 
Haule K, Shim J H and KotUar G 2008 Phys. Rev. Lett. 100 226402 
Miyake T and Aryasetiawan F 2008 Phys. Rev. B 77 085122 

Georges A, KotUar G, Krauth W and Rozenberg M 1996 J, Rev. Mod. Phys. 68 13 
Craco L, Laad M S, Leoni S and Rosner H, arXiv: 0805.3636. 

Shorikov A O, Korotin M A, Streltsov S V, Korotin D M and Anisimov V I, arXiv: 0804.3283; 

JETP 107 N6 (2008), to be pubUshed. 
Kurmaev E Z, Wilks R, Moewes A, Skorikov N A, Izyumov Yu A, Finkelstein L D, Li R H and 

Chen X H, arXiv: 0805.0668. 
KroU T, Bonhommeau S and Kachel T et al. arXiv: 0806.2625. 
Malaeb W, Yoshida T and Kataoka T et al. J. Phys. Soc. Japan 77 093714 (2008). 
Hybertsen M S, Schliiter M and Christensen N E 1989 Phys. Rev. B 39 9028; Hybertsen M S, 

Stechel E B, Schluter M and Jennison D R 1990 Phys. Rev. B 41 11068; McMahan A K, Annett 

J F and Martin R M 1990 Phys. Rev. B 42 6268 
Eskes H and Sawatzky G A 1991 Phys. Rev. B 44 9656 

Maier Th, Jarrell M, Pruschke Th and Keller J 2000 Phys. Rev. Lett. 85 1524; Macridin A, Jarrell 

M, Maier Th and Sawatzky G A 2005 Phys. Rev. B 71 134527; Yin W G and Ku W 2008 J. 

Phys.: Conf. Series 108 012032 
Bocquet A E, Mizokawa T, Morikawa K, Fujimori A, Barman S R, Maiti K, Sarma D D, Tokura 

Y and Onoda M 1996 Phys. Rev. B 53 1161 
Nekrasov I A, Keller G, Kondakov D E, Kozhevnikov A V, Pruschke Th, Held K, VoUhardt D and 

Anisimov V I 2005 Phys. Rev. B 72 155106 
Solovyev I, Hamada N and Terakura K 1996 Phys. Rev. B 53 7158 

Andersen O K 1975 Phys. Rev. B 12 3060 Gunnarsson O, Jepsen O and Andersen O K 1983 Phys. 
Rev. B 27 7144 



Baroni S, de Gironcoli S, Corso A D and Giannozzi P, http://www.pwscf.org 
Wannier G H 1937 Phys. Rev. 52 191 

Marzari N and Vanderbilt D 1997 Phys. Rev. B 56 12847; Ku W, Rosner H, Pickett W E and 
Scalettar R T 2002 Phys. Rev. Lett. 89 167204 
[25] Anisimov V I, Kondakov D E, Kozhevnikov AY et al. 2005 Phys. Rev. B 71 125119 



Coulomb repulsion and correlation strength in LaFeAsO 



13 



[26] Korotin Dm, Kozhevnikov A V, Skornyakov S L, Leonov I, Binggeli N, Anisimov V I and Trimarchi 

G 2008 Europ. Phys. J. B 65 91 
[27] Anisimov V I, Poteryaev A I, Korotin M A, Anokhin A O and Kotliar G 1997 J. Phys.: Condens. 

Matter 9 7359; Lichtenstein A I and Katsnelson M I 1998 Phys. Rev. B 57 6884; Held K, 

Nckrasov I A, Keller G, Eyert V, Blumer N, McMahan A K, Scalettar R T, Pruschke Th, 

Anisimov V I and Vollhardt D 2006 Phys. Stat. Sol. (b) 243 2599 
[28] Hirsch J E and Fye R M 1986 Phys. Rev. Lett. 56 2521 
[29] Vidberg H J and Serene J E 1977 J. Low Temp. Phys. 29 179 
[30] Lie C, Kondo T, Tillman M E et al. arXiv: 0806.2147. 
[31] Korotin M, Fujiwara T and Anisimov V 2000 Phys. Rev. B 62 5696 



